How to optimally parametrize deviations from General Relativity 
in the evolution of cosmological perturbations 
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The next generation of weak lensing surveys will trace the growth of large scale perturbations 
through a sequence of epochs, offering an opportunity to test General Relativity (GR) on cosmo- 
logical scales. We review in detail the parametrization used in MGCAMB to describe the modified 
growth expected in alternative theories of gravity and generalized dark energy models. We highlight 
its advantages and examine several theoretical aspects. In particular, we show that the same set 
of equations can be consistently used on super-horizon and sub-horizon linear scales. We also em- 
phasize the sensitivity of data to scale-dependent features in the growth pattern, and propose using 
Principal Component Analysis to converge on a practical set of parameters which is most likely to 
detect departures from GR. The connection with other parametrizations is also discussed. 



I. INTRODUCTION 

Future weak lensing surveys, like the Dark Energy 
Survey (DES) [1], Large Synoptic Survey Telescope 
(LSST) [2] and Euclid [3] will measure lensing shear and 
galaxy counts at a sequence of redshifts, effectively map- 
ping the evolution of the matter and metric perturba- 
tions. Much like the table top and solar system tests of 
General Relativity [4], this will offer an opportunity to 
verify Einstein's equations that specify the way in which 
the matter, the gravitational potential and the space cur- 
vature are related to each other, and thus test the validity 
of GR on cosmological scales [5] . 

Several parametrizations of modified growth have been 
proposed in the literature [6-14] and they can be sepa- 
rated into two types. The first type [6, 7, 9, 13], which 
can be called "trigger" parameters, are directly derived 
from observations with no need to evolve growth equa- 
tions of motion. They are designed to detect a breakdown 
of the standard model, but their values do not necessar- 
ily have a physical meaning in any theory. The second 
type [8, 10-12, 14] can be called "model" parameters. 
They have physical meanings and unique values in spe- 
cific modified gravity theories and can be used to define 
a consistent set of equations with which to compare the- 
oretical predictions to observations. Theoretical predic- 
tions for these parameters in modified gravity and gen- 
eral dark energy models are studied in [15]. 

Eventually, it will be possible to simultaneously fit a 
given parametrization to a combination of all data that 
probe the growth: CMB, weak lensing, galaxy counts, 
and peculiar velocities. One then needs a system of equa- 
tions that is meaningful across a wide range of scales 
and redshifts. Since trigger parameters are constructed 
out of specific types of observables, they can lead to un- 
necessary complications and inconsistencies if used as 
model parameters for calculating predictions for other 
types of data. Also, the evolution of perturbations on 
super-horizon scales is governed by a set of consistency 
conditions which are separate from the sub-horizon dy- 



namics. Namely, as shown in [17, 18], in the absence of 
entropy perturbations, the space curvature defined on 
hyper-surfaces on uniform matter density, Q, must be 
conserved on scales outside the horizon in order to be 
consistent with the overall expansion of the universe. 
Hence, a consistent system of equations should decou- 
ple the super- and sub-horizon regimes. This separation 
of scales is made explicit in the Parameterized Post- 
Friedmann (PPF) Framework of [11], where a different 
systems of equations are used on super-horizon and sub- 
horizion scales. The advantage of the method advocated 
in this paper and used in MGCAMB [14, 16] is that it em- 
ploys a single system of equations across all linear scales, 
without sacrificing any of the important consistency con- 
ditions. The super-horizon and sub-horizon evolution de- 
couple naturally, without having to be explicitly sepa- 
rated. 

In principle, any modification of gravity can be for- 
mally described in terms of an effective dark energy fluid. 
We show that, by construction, our system of equations 
automatically conserves the energy-momentum of this 
fluid. Furthermore, the same condition that insures con- 
servation of C outside the horizon implies adiabaticity of 
the effective fluid perturbations. 

Aside from the consistency of the modified growth 
parametrization, an equally desired property is its sim- 
plicity. One wants to strike a balance between work- 
ing with the fewest number of parameters possible, yet 
still having enough flexibility in the model to capture 
most of the significant information contained in the data. 
We show how one can determine the optimal number 
of parameters based on a Principal Component Analysis 
(PCA), such as one performed in [19]. 

The rest of the paper is organized as follows. In Sec. II 
we present a system of equations to evolve linear pertur- 
bations in a general theory of gravity, demonstrate their 
consistency across all linear scales, and discuss their in- 
terpretation in terms of an effective dark energy. Sec. Ill 
discusses a strategy for finding an optimal number of pa- 
rameters. In Sec. IV we discuss the connection of our 
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method to other parametrizations in the hterature. Wc 
summarize in Sec. V. 



II. PARAMETERIZED EVOLUTION OF 
LINEAR PERTURBATIONS IN A GENERAL 
THEORY OF GRAVITY 

In order to test gravity against the growth of cosmo- 
logical perturbations in a model-independent way, one 
needs a set of equations to evolve linear perturbations 
without assuming GR. The concept of the universe be- 
ing described by a Friedmann-Robertson- Walker (FRW) 
background metric with small perturbations on large 
scales, and with the matter content that is covariantly 
conserved, is more general than GR and we adopt it as 
our starting point. We consider scalar metric perturba- 
tions about a FRW background in conformal Newtonian 
gauge, for which the line element reads 



-a^ir) [(1 + 2*) dr^ - (1 - 2$) dx^] 



(1) 



where $ and \E' are functions of time and space. We will 
work in Fourier space and, for simplicity, only consider 
cold dark matter (CDM) in our equations. Our discus- 
sion can be generalized to include baryonic and radi- 
ation effects, which would be important at sufficiently 
early times. We assume adiabatic initial conditions and 
covariant conservation of the energy-momentum tensor 
of matter. At linear order, the conservation equations in 
Newtonian gauge are 



3$' = 0, 



(2) 
(3) 



where S is the energy density contrast, v the irrotational 
component of the peculiar velocity, and primes indicate 
derivatives w.r.t. In a. In what follows we will work with 
the gauge-invariant comoving density contrast 
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(4) 



which is particularly convenient on super-horizon scales 
to avoid gauge artifacts. 

In order to solve for the evolution of the four scalar 
perturbations { A, v, $, 4'} we need two additional equa- 
tions, normally provided by a theory of gravity (such as 
GR) which specifies how the metric perturbations relate 
to each other, and how they are sourced by perturbations 
in the energy-momentum tensor. One can parametrize 
these relations as 



where Mp = I/SttG, and ri{a, ky and ii{a, k) are generic 
functions of time and scale; they will assume an ex- 
plicit form once a theory is specified. For instance, equa- 
tions for the growth of perturbations have been derived 
for Chameleon type scalar-tensor models [21] (such as 
f{R) [22, 23] theories) in [24-29], and for the Dvali- 
Gabadadzc-Porrati model (DGP) [30] and its higher- 
dimensional extensions [31] in [32-36]. Given the linear 
perturbation equations for a specific model, it is straight- 
forward to the determine functions 77 and /x [12, 14, 36]. 
In principle, one could also use the functions /x and i] to 
represent dark energy perturbations [37-40] or cosmolog- 
ical effects of massive neutrinos [41]. 

It should be noted that in a general theory of grav- 
ity the superposition principle may not hold. Hence the 
dynamics of large scale density fiuctuations that one de- 
rives by perturbing the background solution at linear or- 
der need not be the same as the result of averaging over 
small scale fiuctuations. A'^-body simulations in scalar- 
tensor theories of chameleon type [43, 43-45], such as 
f(R), and in higher-dimensional extensions of DGP [46- 
48] demonstrated that in those models on large scales 
one recovers the predictions of linear perturbation the- 
ory. However, this need not hold in all modified gravity 
model. With this in mind, one can view /i as a purely 
phenomenological function which relates A and 4* and 
which could be non-linear in a particular modified grav- 
ity theory. In other words, in a non-linear theory fj, may 
itself be a function of A. 

Together, Eqs. (2)-(6) provide a complete system for 
the variables {A,'!;,$,4'}. They can be combined into 
a closed system for the variables A and v. It is in- 
structive to write this system of equations in terms of 
dimcnsionless quantities; specifically, wc will use p = 
SMpHqUmo.'^ and introduce the new variables: 



k 



u — pv 



Em — 



E 



IP 



(7) 



Now, combining equations (2)-(6) and substituting the 
dimcnsionless variables (7) we obtain the following sys- 
tem of equations: 



A' 



2E '/^ 





ivfj-Y 


A + 








H P 



u = — 



IT 



p2 + ^r|^i 

3 E,n . 

u ^AtA . 

2 E ^ 



(8) 
(9) 



Given the functions /i and 77, Eqs. (8) and (9) can be 
integrated numerically to find A and u. The potentials $ 
and 5* can then be determined from (5) and (6). 



77(0, fc), 

„2 



2M|, 



^i{a,k)pA , 



(5) 
(6) 



^ Note that we are now using rj instead of 7 for the ratio of the po- 
tentials in order to avoid confusion with the growth index param- 
eter 7 [6, 7], and to be consistent with other hterature [10, 15, 20]. 
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Eqs. (8) and (9) hold for any k, and can be used to 
evolve cosmological perturbations from super-horizon to 
sub-horizon linear scales if fi{a, k) and 77(0, k) arc pro- 
vided. As we will show, these equations automatically 
ensure that the evolution of perturbations outside the 
horizon is independent of /i up to 0{p'^ j ^v) terms, as 
required by the consistency with the FRW background 
(15). 



A. The ACDM limit 

In the ACDM model, ^ = 1 = rj, and Eqs. (8) and (9) 
simplify to 



A' 



,H' 
'TT 



^ -IE 



3 E,n . 

u A . 

2 E 



(10) 
(11) 



At epochs when the radiation component can be ignored, 
3{H'/H) = -9E,n/i2E) and Eq. (10) becomes 



A' 



(12) 



Combined with (11), this gives the usual second order 
equation 



A" 



IT 



A' 



3 Em 

2ir 







(13) 



which is scale-independent. In a modified gravity, how- 
ever, the time evolution of A will be scale-dependent for a 
general ^ and 77. In particular, it will be scale-dependent 
even if yu and 77 depend only on time, as long as 7^ 0, 
as can be noticed from (8). 



B. Super-horizon evolution 

The consistency of the long wavelength perturbations 
with the FRW background requires that 



C = {5- 3$)' = 0(p2 



(14) 



on super-horizon scales [p — k/{aH) <C 1) for adiabatic 
perturbations [17, 18], where C is the curvature perturba- 
tion on hypersurfaces of uniform density [49] . As pointed 
out in [17], Eq. (14) follows directly from matter conser- 
vation (Eq. (2)), as long as the kv/{aH) term is 0{p^) 
(which is the case in GR and in all viable gravity models 
considered in the literature). The super- horizon conser- 
vation of C implies a second order differential equation 
for the metric potentials [18, 50]: 



TTll 

H' 



H' H" 
ll ^ IT' 



* = 0{p^ 



(15) 



Therefore, once a relation between the two potentials is 
specified, i.e. once is given, Eq. (15) is sufficient to solve 



for the evolution of super-horizon scale metric perturba- 
tions. Then Eq. (14) and the super-horizon limit of Eq. 
(3) can be used to infer 5 and v. 

In a multi-fiuid system, the curvature perturbation on 
the uniforni-total-density hypersurface, Cj remains con- 
stant on super-horizon scales if the uniform-density hy- 
persurfaces for the different fluids coincide on superhori- 
zon scales, i.e. if we have adiabatic initial conditions for 
the multi-fluid system. Since we are mostly interested in 
structure formation during matter domination and later, 
we take ^ = C,„, and analyze the SH behavior of C, in 
our framework (2)- (6). It is straightforward however to 
generalize the following discussion to the multi-fluid case. 

It is easy to see that Eq. (15) follows from the set of 
equations (8) and (9) as long p^ /{^rf) — > in the small p 
limit. Ideed, for p ^ 1, Eqs. (8) and (9) become 



A' = -A 



1-77 (7//i)' 



o H' 
2 + — 
H 



2 H'H 



3 Tj^iEr, 



3 E,n . 

u iiiX 

2 E ^ 



(16) 
(17) 



Combining them into a second order equation for A, and 
using Eq. (6), we obtain 



El i 

IT'^Vj 
IT ^ IF 





H" 


'U 










1 








V. 









(18) 



This equation is equivalent to Eq. (15), after Eq. (5) is 
used to express $ in terms of \E', as long as firj does not 
approach faster than p'^. This requirement is likely to be 
satisfied in any reasonable model, given that fi — r] = 1 
in GR, and radical deviations from GR are typically 
discouraged by the data. Eqs. (2)- (6) are implemented 
in MGCAMB in synchronous gauge, and their super- 
horizon consistency was demonstrated in [14]. 

Eq. (18) shows that in our framework the super- 
horizon evolution of $ and 4* is independent of /i. Eq. (14) 
implies that this is also the case for the evolution of S in 
the Newtonian gauge. On the other hand, ^ and fj,' ap- 
pear explicitly in Eqs. (16) and (17), and it may seem 
that fi could affect the evolution of super-horizon sized 
perturbations. For instance, Eq. (16) implies that if one 
sets 77 = 1 and /i' > 0, the growth of A will be suppressed. 
One way to see that this does not amount to an inconsis- 
tency with Eqs. (15) and (18) is to observe that A itself 
is 0{p^) on super-horizon scales^. Hence, whatever im- 
pact fi has on the super-horizon evolution of A can only 
be of order 0{p^) and, therefore, is unobservable - it is 
completely hidden by cosmic variance. We demonstrate 
this with an explicit example in a separate subsection. 



^ It is well-known that in synchronous gauge 5 is 0{p^) on large 
scales [51]. There the remaining gauge freedom is used to sot 
V = for CDM, so the synchronous gauge 5 is the same as A. 



4 



'=3 0.6 



0,4 



k=0.0001 h/Mpc 

- k=0.001 h/Mpc 

k=0.01 h/Mpc 

k=0.1 h/Mpc 

LCDM 




0.1 



FIG. 1: The growth factor, {A(k, a)/a)/(A(k, ai)/ai), for the 
wavenumbers: fc = 0.1 (dotted hne), 0.01 (long-dashed), 10^'^ 
(short-dashed) and 10~* (dot-dashed) h/Mpc as a function 
of a for ^J,{z) given by (22). The red soUd line is the ACDM 
solution, identical for each k. Note the approximately scale- 
independent enhancement for sub-horizon modes (fc = 0.1 
and 0.01 h/Mpc) at z < 1 due to a rescaling of Newton's 
constant by /x. The long wavelength modes {k — 10^'' and 
k = 10~* h/Mpc) experience a suppression as expected from 
(16) for fj,' > 0. However, because A is 0{p^) for small k, this 
suppression is concealed by cosmic variance as illustrated in 
Fig. 2. 



Sub-horizon evolution 



For p ^ 1 equations (8) and (9) become 

A' = -u 



u = — 



2 H 

H 



3 Em . 

u uA 

2 E ^ 



(19) 
(20) 



They can be combined into a second order equation: 



A" 



ir_ 
IT 



A' 



3 i?„ 
2~B 



■/iA = , 



(21) 



which has the same form as (13) except for the rescaling 
of the Newton's constant by fi. Therefore, the growth of 
A on sub-horizon scales is directly affected by fx(a,k), 
and is independent of 77. 



D. An example 

Let us consider a toy example to illustrate the theo- 
retical arguments made in the previous subsections. We 
choose 77 = 1 and a scale independent form of /i: 



^(z) 



1 



Mo 



1 - flQ 



tanh 



Sz 



(22) 



which describes a transition from fi = 1 at z > Zs to 
fi ~ liQ at z < Zg, with a width set by Sz. We set fiQ ~ 2, 
Zs = I and Sz = 0.1. Note that the transition in fi oc- 
curs on all scales, including those outside the horizon. If 
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FIG. 2: The matter power spectrum at z = for ACDM (black 
solid curve) along with the associated cosmic variance (blue 
shaded region) , and for the modified gravity (MG) example of 
Sec. II D (red dashed curve). Note the potentially observable 
enhancement on sub-horizon scales, while the suppression due 
to the super-horizon variation of jj, is completely hidden in 



cosmic variance. 



we interpret /i as a rescaling of Newton's constant, such 
a super-horizon variation is technically inconsistent with 
the Friedmann equation which sets the relation between 
the background metric and the background density. We 
choose this unphysical example on purpose, to demon- 
strate that, because Eqs. (2)-(6) are by design consis- 
tent with the background expansion on SH scales, the 
super-horizon variation in fj. is undetectable (as expected 
from the arguments in the previous subsections). On sub- 
horizon scales, however, increasing the value of /i should 
result in enhanced clustering. 

In Fig. 1 we plot the gauge-invariant growth factor 



G(fc,a) 



A(fc,a)/a 
A{k,ai)/ai 



(23) 



where a.i is the scale-factor at some initial time. We con- 
sider four values of fc, ranging from 0.1 to 10^** h/Mpc. 
The ACDM solution is scale-independent, as expected 
from Eqs. (11) and (12), and is shown in the red solid 
line. However, the evolution becomes scale-dependent for 
/i(z) ^ 1, and the dashed and dotted blue lines show the 
solution for the four different scales in the case of fi given 
by (22). One can notice an almost scale-independent en- 
hancement of growth for sub-horizon modes (fc = 0.1 and 
0.01 h/Mpc) at z < Zs as expected from Eq. (21). On 
other hand, the modes that were fully (k = lO^^h/Mpc) 
or partially (fc = 10~'^h/Mpc) outside the horizon when 
M began growing, experience a suppression, as expected 
from (16) for /i' > 0. This suppression is not observable 
because A is Oij?^ for small -p and is completely domi- 
nated by cosmic variance. This is demonstrated explicitly 
in Fig. 2, where we plot the gauge- invariant matter power 
spectrum P(fc) oc A^, along with the associated cosmic 
variance. 

On super-horizon scales, the metric perturbations ^ 
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FIG. 3: The evolution of ^{k, a)/'^[k, ai) for for the wavenum- 
bers: fc = 0.1 (dotted hne), 0.01 (long-dashed), 10"^ (short- 
dashed) and 10~* (dot-dashed) h/Mpc as a function of a for 
li{z) given by (22), (as in Fig. 1). The red solid line is the 
ACDM solution, which is scale-independent. Modifying ^ af- 
fects only the modes that cross the horizon. The super-horizon 
gravitational potential, as expected, does not depend on the 
Poisson equation and, hence, the choice of /x. 



and $ determine observablcs such as the CMB tem- 
perature fluctuations. Fig. 3 shows the evolution of 
^(fc, a) /'^{k, ai) for the model (22) studied in this subsec- 
tion. The red solid line is the scale-independent ACDM 
solution. The effect of the transition in /i is profound on 
sub-horizon scales, but becomes smaller on larger scales, 
and eventually disappears for super-horizon scales. The 
independence of the SH evolution of the potentials on /i is 
physically expected; the Poisson equation indeed should 
not play a role on SH scales, where the evolution of ^' 
and $ is uniquely determined as a function of ry (15). 



E. The effective dark fluid interpretation 

In an alternative gravity model, the Einstein equation 
can in general be written in terms of a functional of the 
metric, i.e. E[gfj^,y] = ^nGT^^. In the GR case E[g^y\ 
corresponds to the Einstein tensor G^^. One can always 
interpret the additional terms in E[gf^^] as contributions 
from the energy-momentum tensor of an effective dark 
fluid. That is, formally we can always write 



(24) 



In general, such a fluid will be imperfect, and its scalar 
perturbations will be characterized by a density contrast, 
velocity potential, pressure and a shear stress. Since the 
Bianchi identity is a geometrical property of the Riemann 
tensor, the contracted Bianchi identity 



to associate the contracted Bianchi identity with the dif- 
feomorphism invariance of the action and with the con- 
servation of energy-momentum of matter fields. Indeed, 
applying the contracted Bianchi identity to the Einstein 
equations in GR, one automatically obtains the covariant 
conservation of the energy-momentum tensor. Namely, 
G^^ = SttGT^^ together with (25) implies £>T^^ = 0. 
Similarly, if we assume D^^T^^ — O'^, as we did earlier, 
then the effective fluid will also be covariantly conserved 
at all orders in perturbation theory. 

With two equations provided by the conservation of 
T^[J, at the linear level in scalar perturbations we have 
two remaining degrees of freedom required to specify the 
evolution of the dark fluid. These can be interpreted 
as the shear stress and the sound speed, which in our 
parametrized formalism of Eqs. (2)-(6) can be determined 
once the functions /i(a, k) and ry(a, k) are specified. By 
construction, will be conserved for any choice of /i 
and 77 and the evolution of perturbations in the effective 
fluid will be completely determined. 

While the conservation of does not impose any 
conditions on /z and 1], we next show that there will be 
a constraint if we further require perturbations of the 
effective fluid to be adiabatic. Using the version of (24) 
linear in scalar perturbations, the effective energy density 
and pressure perturbations can be written as 



2Af2fc2 

t: $ 



6Af|i?2 ($' + , (26) 



SP, 



off 



-6P + 2MlH'^ 



2^ + 3 1 * 



IF 



3 $' + 



($ - *) 



(27) 



where p5 and 6P are respectively the density and pressure 
perturbation of matter flelds. Using Eqs. (5) and (6) we 
could further express these effective quantities in terms of 
the matter variables, 77 and /x, and their evolution would 
be fully specifled for given 77 and /.t. 

Taking the super- horizon limit (p <C 1) wc flnd 



5P 



off 



-a(off) 



$"4 



TTtI 

H' 



El 
IF 



(28) 



where c^^^^g-^ is the adiabatic speed of sound of the effec- 
tive fluid, i.e. 



''a(cff) 



Wcff 



''off 



3(1 + Woff) 



(29) 



-D^G^, ^ D^' [R^, - i<?^.i?j = (25) 

holds independently of the form of the gravitational ac- 
tion at all orders in perturbation theory. It is common 



^ This assumption in part follows from our choice to work in the so- 
called Jordan frame, in which matter fields follow the geodesies 
of the metric. In other words, it amounts to the assumption that 
it is possible to write an action in which the matter Lagrangian 
is minimally coupled to the metric. 
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FIG. 4: The three best constrained eigenmodes of /i for DES, 
adapted from Ref. [19] 



and where we have assumed that matter and the effective 
fluid have adiabatic modes on long wavelengths, i.e. that 
they share a uniform-density hypersurface. 

From Eq. (28) we see that in order for the effec- 
tive fluid to have adiabatic pressure perturbations, i.e. 



SPoS 

have 



— ^ 



c2icff)^PcS on super-horizon scales, one must 



El 
IF 



EL 
IT 



El 
IT' 



* = , (30) 



which is equivalent to the consistency condition (15), de- 
rived from the conservation of the curvature perturbation 
C, in absence of entropy perturbations. In the previous 
subsection we have demonstrated that Eqs. (2)- (6) sat- 
isfy the super-horizon conservation of C as long as 
goes to zero in the p — >■ limt. 



III. HOW MANY PARAMETERS TO FIT? 

In the previous section we showed how two functions, 
fj, and T], can be used in a consistent way to parametrize 
the linear growth of perturbations in a general modifi- 
cation of gravity. We have not, however, discussed the 
actual parametrization of the functions themselves. This 
is not an issue when testing a particular theory, where 
fi and 77 have a specific time- and scale-dependence and 
can be specified with just a few parameters, e.g. via a 
direct use of Eqs. (36) or their equivalent [12, 14]. How- 
ever, one may want to measure fj. and 77 from the data 
without necessarily assuming a particular class of mod- 
els. The question then is how to strike a balance between 
simplicity, i.e working with as few parameters as possi- 
ble, and allowing for enough fiexibility in these functions 
to capture all of the information contained in the data ^. 

For example, one can discretize ^(fc, z) and r]{k, z) on 
a grid in {k, z) space and treat the values of the func- 



* The discussion of this section is specific to model parameters and 
does not apply to triggers, which will be considered in Sec. IV. 
Triggers are reconstructed directly from data, so there is no con- 
trol on their form. However depending on their definition they 
might capture or not some important features in the data. 




FIG. 5: The best measured eigenmodes of w from SNe (black 
solid), CMB (red dash), galaxy counts (green dash dot), weak 
lensing(blue dash dot dot) and combined (magenta short 
dash). This forecast result is adapted from Ref. [54]. 



tions at each grid point, which we will call pixel, as inde- 
pendent parameters. In some of the earlier literature, /i 
and 77 (or a related set of parameters) were taken to be 
constants, or have two very wide pixels parameterizing 
a transition form GR at early times, to different values 
today [20, 52]. Scale-dependent variations have typically 
not been considered. This was motivated primarily by 
simplicity. Namely, the idea is to start with the simplest 
possible model and complicate it only if the fitted values 
of this simple model parameters shows hints of departure 
from ACDM. However, this logic may not always be ap- 
propriate in modified gravity studies and might lead to 
missing on important information contained in the data. 

Indeed, as was clearly shown in [19], the growth data 
is much more sensitive to the shape of the functions 
/i(fc, z) and ry(fc, z), especially to their fc-dependence, than 
to their overall amplitude (or time-dependence). In [19], 
a Principal Component Analysis (PCA) was performed 
to find the best constrained uncorrelated linear combi- 
nations of the ^(fc, z) and ri{k, z) pixels ~ the so-called 
eigenmodes ~ for several surveys. For example, the three 
eigenmodes of best constrained by the Dark Energy 
Survey (DES) [1] are shown in Fig. 4. Note that the best 
constrained mode has a node along the fc-direction, i.e. it 
is definitely not well approximated by a constant function 
of scale. More generally, in [19] the authors have found 
that the well-constrained modes of fj. and rj for several 
surveys display a non trivial pattern in the k direction, 
with no eigenmodes that have no nodes. 

It is instructive to compare this with the findings of 
a PCA analysis on the dark energy equation of state 
w{z) [53; 54]. In general, w is a time-dependent function, 
and the simplest parametrization is that of a constant w. 
As it was found in [53, 54], almost every best constrained 
mode of w forecasted for different surveys has no nodes 
in z, as shown in Fig. 5. This means that observables are 
indeed most sensitive to an average value of w and one 
can expect the tightest constraints on w{z) when fitting 
a constant w, with constraints on variations in w being 
weaker. 

This is not the case for the functions fi{k, z) and 
r]{k, z), especially in regard to their scale-dependence. Be- 
cause of degeneracies with other parameters such as the 
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rJm, w, and as, the growth observables are rather insen- 
sitive to the average values of z) and r]{k, z) over 
scale, and not too sensitive to their time-dependence. 
On the contrary, data is rather sensitive to their scale- 
dependence, as shown clearly by the modes of Fig. 4. This 
is not surprising - for example, redshift surveys primar- 
ily constrain the shape of the power spectrum P{k), with 
the amplitude being degenerate with other parameters. 

The above arguments prove that when fitting /x(fc, z) 
and r](k, z) to data, it is important for their parametriza- 
tion to be flexible enough to allow for some scale- 
dependence. Otherwise, we would not be exploiting the 
true discovery potential of the data, and might risk to 
miss on detecting departures from ACDM. 

Still, while allowing for scale-dependence, one needs to 
choose an optimal pixelization, and the PCA method can 
be of help once again. Indeed, it allows to determine the 
optimal number of pixels necessary to extract most of the 
information from the data. Having too many pixels makes 
the fitting to data numerically difficult, while having too 
few may lead to misleading results. A reasonable way 
to arrive at the optimal number, would be to perform a 
Fisher Matrix forecast using the specifics of the particular 
survey, such as DES and LSST, and then perform a PCA 
on the Fisher matrix. The best constrained eigenmodes 
can then suggest the minimal number of pixels needed 
to describe them. For example, the eigenmodes shown in 
Fig. 4 suggest that one should allow for the pixelization of 
fi to have at least two fc-intervals, so that it can describe 
a potential variation at ^0.1 h~^Mpc. They also suggest 
that data cannot constrain sharp transition in time, and 
one can have a single redshift pixel, with the transition 
at z ~ 1.5. 



IV. CONNECTION WITH OTHER 
PARAMETRIZATIONS 

Several alternative parametrizations of modified linear 
growth have been proposed in the literature and they 
can be divided into two types. The first type consist 
in parameters that can be called triggers [6, 7, 9, 13]. 
They can be derived directly from observations, with no 
need to evolve growth equation. Any disagreement be- 
tween the observed trigger parameter and its value in the 
ACDM model would indicate some sort of a modification 
of growth. The second type consists of parameters and/or 
functions which can be called model parameters and arc 
introduced to play a role similar to our functions fi and rj. 
Namely, they are used to build a system of equations for 
the evolution of metric and density perturbations. Un- 
like the trigger types, which are directly calculated from 
the observables, the model parameters are measured by 
fitting the model that they define to data. Both types of 
parameters have their purpose. Working with functions 
like ^ and rj gives one a consistent set of equations with 
which to compare theoretical predictions to observations. 
A measurement of a particular form of scale- or time- 



dependence of these functions could directly rule out, or 
point to a theory. The trigger parameters, on the other 
hand, are designed to detect a breakdown of the stan- 
dard model, but their value does not necessarily have a 
physical meaning in any theory. 

For example, the commonly used trigger parameter 
7 [6, 7] is defined via 



din a 



In 



A(fc,a) 
A{k,a,) 



(31) 



where ^^(a) = Em/E. As shown in [6], 7 — 6/11 w 0.55 
provides a solution to Eq. (13) if one neglects terms 
C[(l-/)^] (orO[l-0„(a)]2). Since most of the informa- 
tion on linear clustering is expected to come from epochs 
before the matter-A equality, the approximation involved 
in this parametrization is not unreasonable. In principle, 
/ can be extracted by observing clustering at several red- 
shifts, while 51m (a) can be measured from background 
expansion probes, like CMB and supernovae. Then a sig- 
nificant deviation of the observed 7 from its predicted 
value of 0.55 would indicate a breakdown of ACDM. 

The trigger parameters, such as 7 can in principle be 
used also as model parameters. For example, given 7, one 
can integrate Eq. (31) to find A. To find the gravitational 
potentials $ and their relation to each other and to A 
must be specified additionally. In other words, one needs 
to provide a total of three functions instead of two (e.g. 
fi and T] introduced earlier), or assume that one or more 
of the Einstein equations are valid. By working with 7, 
one essentially no longer assumes conservation of matter 
energy-momentum. Alternatively, one can use Eqs. (21) 
(since 7 is only meant to be used to characterize the 
growth on sub-horizon scales) and (31) to express fi in 
terms of 7: 



+ 77 



H 



7' ln(r2m) 



(32) 



For a flat ACDM background eq.(32) simplifles to 



M = -fir' [^7a + '2- 37 



+ 3(7-i)ll„+7'ln(r!„)] . 

(33) 

Then, for a given 7, and with an additional speciflca- 
tion of 77, one can use (32) in Eqs. (2)- (6) to flnd con- 
sistent solutions for the linear perturbations. Note that 
with 7 = 6/11 substituted into (33), fj, is not a constant, 
but a slowly varying function evolving from fi = I dur- 
ing matter domination to /i ~ 1.04 today (for standard 
ACDM parameters) due to the 0[{1 — /)^] corrections. 

Given that trigger parameters are designed in terms 
of a specific observable and for a limited purpose, they 
can lead to unnecessary complications and inconsisten- 
cies if used as model parameters for calculating predic- 
tions for other types of data. When performing a global 
fit to all available data, one needs a consistent system 
of equations for calculating predictions for all of the ob- 
servables: CMB, weak Icnsing, galaxy counts and peculiar 
velocities. 
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As discussed above, in addition to the conservation 
of energy-momentum one needs to specify two functions 
relating A to the potentials, and the two potentials to 
each other. The actual choice of the two relations that 
define such two functions is no unique. However, work- 
ing with ?7 (5) and /i (6) has certain advantages. One is 
that on super-horizon scales only ?/ affects the pertur- 
bations, while fi naturally becomes irrelevant. Also, on 
sub-horizon scales, only fi affects the growth of matter 
over-densities. The physics determining the evolution of 
perturbations on super-horizon scales is not necessarily 
related to their sub-horizon dynamics. This distinction 
is made explicit in the Parameterized Post-Friedmann 
(PPF) Framework of [11], where different systems of 
equations are used on super-horizon and sub-horizon 
scales. On super-horizon scales one must provide only 
g = (1 — ?/)/(l + r/), while on sub- horizon scales one needs 
g and /g = -1 - a2pA/(M|,fc^($ -I- <i>)). A transition 
between the two regimes needs to be additionally speci- 
fied in a way that is consistent with the conservation of 
energy-momentum. The advantage of the PPF is that it 
does not need even weak assumptions on g and fa in or- 
der to satisfy the consistency conditions, while our fi and 
r] technically must obey p^/^^rj) — >■ for small p. The 
latter, however, is an extremely mild assumption, that is 
likely to be satisfied in any reasonable modified gravity 
theory. In turn, the benefit of our approach is that one 
evolves the same system of equations on all linear scales, 
with no need for an additional transition function. 

Sometimes, it is convenient to define the function 
I](fc,a) 



S(a, k) 



(1 + ^) 



(34) 



which is directly related to the lensing potential $ -f- ^ 
(much like fa in PPF). As such, weak lensing (WL) mea- 
surements are sensitive to S as well as the Integrated 
Sachs- Wolfe effect (ISW), which is determined by the 
time variation of the lensing potential. Thus E is more 
directly constrained by these observations than or ^. 
Another advantage of this function is that in the popular 
scalar-tensor/Chameleon theories, as well as in higher- 
dimensional theories, such as DGP [30], this function as- 
sumes a simple expression and is effectively 1 on small 
scales [15]. Specifically, in Chameleon type theories in 
which the scalar field couples to CDM with a coupling 
a(0) one has 



fj,{a,k) 



r]{a,k) 



ni^a^ + (l 



W) k^ 



m^a^ + (1 - \oi^^) e 
ni^a^ + (1 + 



-q(0)/Mf 



(35) 



(36) 



where = da/dc/j and m is a time-dependent effective 
mass scale. Note that both fi and 77 are functions of time 
and scale. On the other hand, E is only a function of 
time: 



S(a, k) 



In the DGP model and its higher dimensional generaliza- 
tions we have [36] 



S(a, k) 



1 



1 + {a/kr. 



(38) 



where is a characteristic lengthscale of the model, while 
a = 1/2 for DGP and for two or more extra-dimensions. 
Since on super-horizon scales one function is sufficient to 
fully describe the evolution of perturbations, it is best to 
use S in pair with rj^ and not with fi. This way, it will 
automatically becomes negligible on super-horizon scales 
and 77 will be the only important function. 

In Caldwell et al. [8] the authors introduced what they 
called the gravitational slip, i.e. a function zu parametriz- 
ing the difference between the gravitational potentials as 



f 1 - 

■CD — 1 = 



V 



(39) 



It should be used in combination with another function 
parametrizing the relation of metric potentials to matter 
density contrast, e.g. fi. Also, generally one should allow 
for it to be scale-dependent. 

Amendola et al. [10] use a function equivalent to zu 
(39) in combination with either E (34) or Q, defined as 



fc2A/2$ 2E?7 
— -Yir- = l"l = 



pa^A 



1 



ri 



(40) 



(37) 



As for E, an appropriate pair to describe modified growth 
would be (Q, r]). 

Different observables will probe different combinations 
of the density contrast and metric potentials correspond- 
ing to one or more of the parametrizing functions. For 
example clustering, i.e. the growth of A, responds to the 
potential ^, and therefore to the function /i. However, it 
is also affected by the magnification bias, i.e. the lensing 
of the clustering pattern of the background sources by 
the intervening gravitational potentials. Hence, any real- 
istic observation of clustering will be sensitive also to E 
(and hence ry) as well as pL. Weak lensing and the ISW 
effect probe the sum of the potentials ($-1-^) and, there- 
fore, are directly related to the function E; of course they 
also constrain p and 77, with stringent constraints on pi as 
shown in [19]. Finally, peculiar velocities respond directly 
to the potential 4', therefore to the function /x. 

In summary, given that two functions are necessary to 
parametrize departures from LCDM in the growth pat- 
tern, one needs to decide which pair among the several 
functions described above, is best suited for a given anal- 
ysis and for presenting the corresponding results. Our 
choice of the pair (77, pi) has been guided by the fact that 
on super-horizon scales p naturally becomes irrelevant 
and we are left with only one function, -q, as expected on 
those scales. Notice also the importance of p being de- 
fined through the Poisson equation involving A and not 
(5; this allows for p to be defined consistently on all scales 
and to go exactly to unity on all scales for LCDM. From 
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[r], fi), we can derive other parameters which may be 
more suitable for interpreting observational constraints. 
For example, wc can choose to describe departures from 
GR in terms of (E, fi) because WL measurements and 
the ISW effect are sensitive to S, and fi can be de- 
termined from peculiar velocity measurements, while rj 
is not probed directly by any observable and would be 
highly degenerate with /x or S. 

V. SUMMARY 

As of today, the cosmological concordance model, 
ACDM, provides a good fit to all available observations. 
However, upcoming and future weak Icnsing surveys will 
bring a significant improvement in cosmological data sets 
and will offer an unprecedented opportunity to test GR 
on cosmological scales. As demonstrated in [19], the cur- 
rently available clustering information pales in compar- 
ison to what we will learn from a survey like LSST. A 
number of alternative gravity theories that are currently 
indistinguishable from ACDM will be tested, potentially 
providing clues about the physics causing cosmic acceler- 
ation. Even if no hints of new physics are observed, at the 
very least we can directly confirm the validity of GR at 
epochs and scales on which it has not been tested before. 

This paper is a step towards building an optimal frame- 
work for testing GR. We showed how a single system of 
equations can be used consistently to evolve perturba- 
tions across all linear scales. These equations are imple- 
mented in MGCAMB [14, 16], a publicly available mod- 
ification of CAMB [55], which can be used to evaluate 
CMB, weak Icnsing, number counts, and peculiar veloc- 
ities spectra for a choice of functions fi and rj. We have 
also argued, based on the scale- and redshift-dependent 
patterns of the best constrained eigenmodes, that fu- 
ture surveys will primarily probe the scale-dependence, 
and not so much the overall normalization or the time- 
dependence of these functions. Hence, in order to fully ex- 



ploit the discovery potential of data, parametrized mod- 
ifications of gravity must allow for scale-dependence. 

While we have not addressed it in this paper, the 
galaxy counts only trace the underlying density field 
up to a bias factor. As explained in [56], any scale- 
dependence in the linear growth, will result in a scale- 
dependence of the linear bias. The latter however, is 
directly related to the scale-dependence of A{k, z) and 
therefore can be determined once /Li(fc, z) is specified. 

Most of the clustering information comes, and will con- 
tinue to come, from scales that have crossed into the non- 
linear regime, not described by the linear parametrization 
of Sec. II. One could test gravity on non-linear scales 
by designing trigger parameters that would indicate a 
breakdown of GR. Alternatively, A'^-body simulations and 
higher order pcrturbativc expansions [57] can be used to 
constrain particular types of modified gravity theories 
[58]. 

The precise accuracy of the tests will strongly depend 
on the ability to control the systematic effects, and their 
extent will not be fully known until experiments begin 
operating. Some preliminary estimates of the effect of 
systcmatics on cosmological test of GR have been re- 
ported in [19] and a more comprehensive analysis will be 
presented in [59]. 

Note added: while we were preparing the manuscript 
for submission, a paper appeared on arXiv [60] where 
some of the issues raised in this paper are also discussed. 
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